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We have been developing a program package called GRACE/SUSY-loop which automati- 
' cally calculates the MSSM amplitudes in one-loop order. We present numerical results 

(~| , of calculations for pair- production and three-body decay of the lighter stop (ti) at the 

International Linear Collider (ILC) using GRACE/SUSY-loop. Since the distributions 
of missing transverse energy (MET) depend on mass spectrum of SUSY particles, we 
consider two scenarios on three-body decay of ti. In these scenarios, both QCD and 
EW corrections have positive sign for decay widths and cross sections. 



1 Introduction 



. Supersymmetry (SUSY) between bosons and fermions at the unification-energy scale is 

' considered one of the most promising extension of the standard model (SM) of particle 

physics. Among the supersymmetric theories, the minimal supersymmetric extension of the 
SM (MSSM) is a wcU studied framework of SUSY. Since SUSY is a broken symmetry at the 
electroweak-cnergy scale, its relic is expected to remain as a rich spectrum of heavy SUSY 
particles, which makes calculation of amplitudes very complicated. 

For many-body final states, each production process or decay process is described by a 
large number of Feynman diagrams even in tree-level order. There are still more Feynman 
diagrams in one-loop order even for two-body final states. For this reason, we have de- 
veloped the GRACE system [l], which enables us to calculate amplitudes automatically. A 
version of the GRACE system called GRACE/SUSY-loop [2] is a program package for automated 
computations of the MSSM in one-loop order. Other groups independently developed pro- 
gram packages SloopS |5] and FeynArt/Calc [2] for automatic calculations of the MSSM 
amplitudes in one-loop order. 

Recently, we have calculated radiative corrections of production processes and decay 
processes of SUSY particles in the framework of the MSSM using GRACE/SUSY-loop 
In this paper, we show numerical results of the MET distributions for pair-production and 
three-body decay of ti at the ILC in tree-level order, and also show those of the radiative 
corrections to cross sections for pair-production [7] and decay widths for three-body decay 
of ti [Slin] in one-loop order. 
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2 MET distributions 



In the experiments at the Large Hadron CoUider (LHC), searches for SUSY particles are 
performed using an experimental cut on MET as is required to be larger than 100 GeV [TU] 
because of background-events avoidance. When a SUSY particle is heavy and decays directly 
to a charged particle and the lightest SUSY particle (LSP), e.g. — >■ t + Xi, the MET is 
expected to be large because the LSP is not observed in detectors and the momentum of 
the parent particle is small. It, however, can be smaller if a SUSY particle is light and has a 
large momentum because momenta of particles produced from the decay arc boosted along 
that of the parent particle. Although no light SUSY particle has been observed at the LHC, 
there is room for it to survive if MET is small. 

2.1 Two scenarios 

Since the Yukawa coupling of top quark is large, the mass of ti can be smaller than that of 
other squarks. So here we focus on ti pair-production + e+ — > ti + i* and decay at the 
ILC. To investigate the possibility that M^^ is small enough, we set the mass range of ti 
below the threshold of the two-body decay. Since the MET distributions depend on mass 
spectrum of SUSY particles, we consider two scenarios for three-body decay of ti as follows: 

Scenario 1. Large slepton masses 

The major decay mode is ii — ?> bW'^Xi [SlE]j and its branching ratio is nearly 100% if 
the mass range of ti is between the threshold of the three-body decay, Mj_^ ~ 279 GeV 
and the lowest threshold of the two-body decay, which is for the mode ti t + Xi: 
Mj = 368 GeV. 

Scenario 2. Small slepton masses 

'Semi-Leptonic' decay modes, ti -> bl^Di and ti — > bl^vi (l ~ e,ii,T), dominate, and 
the branching ratio of the mode ti — > bW^Xi negligible small due to the small phase 
space. The lowest threshold of the two-body decay, which is for the mode ti — )• b + xt , 
is Mj^ = 298 GeV. 

We set SUSY parameters for these scenarios as in Table [TJ 
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Table 1: SUSY parameters for two scenarios. 
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2.2 Numerical results of MET in tree-level order 

In Scenario 1, BR{ti bW~^Xi) — 1 in the range of for which three-body decay of ti 
is dominant as studied in [5l|6]. In Scenario 2, however, there are several 'Semi-Leptonic' 
decay modes. Figure [T] shows numerical results of Mj_^ dependence of branching ratios for 
three-body decay of ti in tree-level order according to Scenario 2. 
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Figure 1: Mj_^ dependence of branching ratios for three-body decay of ti in tree-level order 
according to Scenario 2. 

Figure [2] shows numerical results of (a) MET distributions for + e+ ti + ^• 
(bW^Xi) + ibW~Xi) according to Scenario 1 and (b) those for e~ -|- e+ — >■ -I- — )■ 
{be'^'ue) + (be~i'e) according to Scenario 2, respectively, at ^/s = 1 TeV in tree-level order. 
In our calculations, events are generated by SPRING [H] which is a built-in package in the 
GRACE system. In Scenario 1, the MET is calculated for the two LSP's. In Scenario 2, a 
major decay mode ti ^ be^ Ve is selected, and the MET is calculated for the two sneutrinos. 




Figure 2: MET distributions at y/s = 1 TeV in tree-level order, (a) for + e+ ti + 1 j ^ 
{bW+Xi) + (bW~Xi) according to Scenario 1; and (b) for e" -I- e+ -> <i -I- (&e+i^e) + 
{be~vl) according to Scenario 2. Red lines indicate Mj^ = 280 GeV. Green lines in (a) and 
(b) indicate Mj^ = 300 GeV and 250 GeV, respectively. 

We should note that the MET in Scenario 1 is small because the difference of numerical 
values between M^^ and the sum of Mw and Af^o is small. Although each W boson decays 
to a pair of quarks or that of a charged lepton and a neutrino, resultant MET distributions 
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Figure 3: Mj-^ dependence of decay widths for three-body decay processes of ti and their 
correction rates according to Scenario 2. 
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are similar to Figure [5] (a) except for the normalization factors because the momenta of the 
final particles from the W-boson decay are boosted along that of ti. It is also notable that 
the peaks of the MET distributions are less than 100 GeV in Scenario 2. 

3 Numerical results of radiative corrections in one-loop order 

In Scenario 1, the radiative corrections to decay widths of ti in one- loop order have already 
been studied in [5||6]. In Scenario 2, decay widths for several 'Semi-Leptonic' decay modes of 
ti are in the same order. Figure [3] shows numerical results of A/^^ dependence of decay widths 
of ti and their correction rates according to Scenario 2. For the decay modes ti bl^vi 
[(a) - (d) in Figure |3], (^Fqcd > (^Few > 0, and for the decay modes ti — >• hl^vi [(e) - (j) 
in Figure [5], (5Few > (^rqcD > 0. The difference in the behavior of 5Few between the two 
is due to that of the QED corrections. Each decay width is enhanced by the total radiative 
correction in one-loop order. 

Figure H] shows numerical results of dependence of branching ratios for three-body 
decay of ti in one- loop order according to Scenario 2. The results are similar to those in 
tree- level order (Figure [1} because (5F's of the decay modes are in the same order. 



Figure 4: Mj_^ dependence of branching ratios for three-body decay of ti in one-loop order 
according to Scenario 2. 

Figure [5] shows numerical results of (a) Mj_^ dependence of cross sections for the ii-pair 
production and succeeding decay process, + e+ ^ <i + t| — > {bW^Xi) + {bW~Xi) 
(be'^h'eXi) + {be~i'eXi) according to Scenario 1 and (b) that for the ti-pair production and 
succeeding decay process, e~ + ^ ti + tl {he'^Ve) + (be~i'e), according to Scenario 2 
at y/s = 1 TeV. For comparison, a decay mode e+ is selected in Scenario 1, and 

a decay mode — >■ fee"*" Ve is selected in Scenario 2. In both scenarios, Jctqcd > Sa-Ew > Oj 
and cross sections are enhanced by the total radiative corrections in one-loop order. 

Wc should note that the MET distributions in one-loop order are similar to those in 
tree-level order except for the normalization factors, and that the MET distributions in the 
two scenarios can be distinguished by the bbe~ tagging. 

4 Summary 

We have calculated the MET distributions for the pair-production and three-body decay of 
ti in tree-level order, and also calculated the radiative corrections in one-loop order of the 
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Figure 5: Mj_^ dependence of cross sections at ^/s ~ 1 TeV. (a) for e + e+ ^ ti + 

tl ^ (bW+Xi) ^{bW~Xi) (^e+i/eX?) + (be'DeXi) according to Scenario 1; and (b) for 
+ ^ ti + tl ^ (fee+i^e) + {he~v'^) according to Scenario 2. 



MSSM for two scenarios on slepton masses. We found that the peak of the MET distributions 
can be less than 100 GeV if M^^ < 300 GeV in both scenarios, so events for MET < 100 
GeV should be also analyzed in detail for ILC study. In the two scenarios, both QCD and 
EW corrections in one-loop order have positive sign for decay widths and cross sections in 
the range of Mj_^ for which three-body decay of ti is dominant. 
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